Spontaneous Relaxation of a Charge Qubit under Electrical Measurement 
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In this work we first derive a generalized conditional master equation for quantum measurement 
by a mesoscopic detector, then study the readout characteristics of qubit measurement where a 
number of new features are found. The work would in particular highlight the qubit spontaneous 
relaxation effect induced by the measurement itself rather than an external thermal bath. 
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The recent renewed interest of measuring a two-state 
quantum system (qubit) stems largely from the rapidly 
developing field of quantum computing. A possible solid- 
state implementation of such measurement is to measure 
a charge qubit by a mesoscopic detector which, for in- 
stance, can be either a quantum-point-contact (QPC) 
Em II 3, or a single-electron-transistor (SET) Q. 

For a realistic setup of such measurement, the non- 
trivial correlation between the detector and the qubit 
has been the focus of recent theoretical studies. How- 
ever, sometimes the treatment of this correlation is in- 
complete. For instance, in a number of publications on 
the qubit measurement by a QPC [ESIEQ) the energy 
transfer between the detector and qubit has been ignored, 
which leads the qubit to an incorrect statistical mixture 
under low measurement voltage, as shown in our work 
Two recent publications considered the energy-exchange 
induced inelastic effect on the detector power spectrum 
by usingj^respectively, the real-time Green's function ap- 
proach Q and the quantum jump technique where 
a number of controversial results were arrived and cause 
further debate plj . 

In thispaper, by generalizing the work of Gurvitz et 
al. 0, 0, 3 , we present an alternative approach to study 
the inelastic effect in the qubit measurement by a QPC. 
Connections with the previous work will be established in 
a transparent way, and new features will be illustrated in 
both the output current and power spectrum. In partic- 
ular, we shall highlight the qubit spontaneous relaxation 
effect induced by energy exchange with the measuring 
device, instead of coupling to an external thermal bath 
as discussed in Ref. 0. 

For the sake of generality, we first formally consider an 
arbitrary quantum system measured by a QPC, described 
by 

H = Ho + H', (la) 
Ho = Hs + Y,{e^clck + e^'dldk), (lb) 

k 

H' = 5][T^9fe{lV's)(^s|}4cfc-t-H.c.]. (Ic) 

k.q 

In this decomposition, the free part of the total Hamil- 
tonian Ho contains the Hamiltonians of the measured 
system Hg and the QPC reservoirs (the last two terms). 



The interaction Hamiltonian H' describes electron tun- 
neling through the QPC, e.g., from state |fc) in the left 
reservoir to state \q) in the right reservoir, with tunnel- 
ing amplitude Tqfc{|'(/'s)('0s|} that is conditioned by the 
eigenstate \ips) of the observable. 

Regarding the tunneling Hamiltonian H' as perturba- 
tion, on the basis of the second-order cummulant expan- 
sion we can derive a formal equation for the reduced den- 
sity matrix as E3 



p(t) = -z/:p(i) - / dT{c'{t)g{t,T)c'{T)g^{t,T))p{t).{2) 



Here the Liouvillian superoperators are defined as 
£(•••) ^ [Hs,{---)], £'(•••) ^ [H'.{---)], and 
g{t,T){---) = G(i,T)(---)Gt(t,r) with G{t,T) the usual 
propagator (Green's function) associated with Hg. The 
reduced density matrix p{t) = Tr£)[pT(i)], resulting from 
tracing out all the detector degrees of freedom from the 
entire density matrix. However, for quantum measure- 
ment where the specific readout information is likely to 
be recorded, the average should be performed over the 
unique class of states of the detector we are trying to 
keep track of. 

For the measurement setup under study, the relevant 
quantity of readout is the transport current i{t) in the 
detector, or equivalently, the number of electrons that 
have tunnelled through the detector, n(t) = J^dt'i{t'). 
We therefore classify the Hilbert space of the detector 
as follows. First, we define the subspace in the ab- 
sence of electron tunneling through the detector as V^^^ , 
which is spanned by the product of all many-particle 
states of the two isolated reservoirs, formally denoted as 
X>(") = span{|*L) (g) l^-ii)}. Then, we introduce the tun- 
neling operator p /Jj, = djc^, and denote the Hilbert 
subspace corresponding to rt-electrons tunnelled from the 
left to the right reservoirs as — (/t)"I?(*^), where 

n — 1,2, ■ ■ ■ . The entire Hilbert space of the detector is 
p = ©„!?("). 

With the above classification of the detector states, 
the average over states in V in Eq. ((2Jl is replaced with 
states in the subspace P^"-*, leading to a conditional mas- 
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ter equation 

JQ 

x£'{r)gHt,T)pT{t)]. (3) 

Here p'^'^\t) = Tr£,(„) [^^(i)], which is the reduced den- 
sity matrix of the measured system conditioned by the 
number of electrons tunnelled through the detector un- 
til time t. Now we transform the Liouvillian operator 
product in Eq. Q into the conventional Hilbert form: 

c'{t)g{t,T)c\T)g\t,T)pT{t) 

= [H'{t)G{t,T)H'{T)G\t,T)pT{t) 

-G{t,T)H'{T)G^{t,T)pT{t)H'{t)]+ll.c. 
= [I - II] + H.c. (4) 

For the convenience of description, we rewrite the in- 
teraction Hamiltonian as II'{t) = QF{t). Here we have 
assumed the tunneling amplitude T^q to be real and inde- 
pendent of the reservoir-state "fcg" , and denoted it by Q 
which depends on the state of the measured system. The 
detector fluctuation is described by F{t) = fit) + f\t), 

with / = Y.kq A'^a ^^"^ = Y.kq ^l^k- To proceed, two 
physical considerations are further involved as follows: (i) 
Instead of the conventional Born approximation for the 
entire density matrix prit) — p{t) ® po, we propose the 

ansatz pT{t) — X^n ® Pd^^ where p^' is the den- 

sity operator of the detector reservoirs with n-electrons 
tunnelled through the detector. With the ansatz of the 
density operator, tracing over the subspace yields 

Tr^(„,[/] - Tro[F{t)Fir)p^S^] 

x[QG(t,r)gGt(t,T)p("'] (5a) 
Tr^(„,[//] = Trn[fHr)p^r'^f{t)] 

x[Git,r)QGHt,T)p'^''-^^Q] 

+TTn[f{r)p^S^'^fHt)] 
x[G(t,r)QGt(i,r)p("+i)g]. (5b) 

Here we have utilized the orthogonality between states 
in different subspaces, which in fact leads to the term 
selection from the entire density operator px- (ii) Due to 
the closed nature of the detector circuit, the extra elec- 
trons tunnelled into the right reservoir will flow back into 
the left reservoir via the external circuit. Also, the rapid 
relaxation processes in the reservoirs will quickly bring 
the reservoirs to the local thermal equilibrium state de- 
termined by the chemical potentials. As a consequence, 
after the procedure (i.e. the state selection) as done in 

Eq. ||SJ), the detector density matrices p^'' and p^^^"* 

in Eq. lO can be well approximated by p^'', i.e., the 
local thermal equilibrium reservoir state. Under this 
consideration, the detector fluctuation correlation func- 
tions become, respectively, {f^{t)f{T)) — C^'^^t — r), 
(/(t)/t(r)) = G(-)(t-r), and {F{t)F{T)) = C{t - 



t) = C^+\t - r) + C^-\t - t). Here, (• • •) stands for 

T:td[{---)p%\ 

Under the Markovian approximation, the time integral 
in Eq. Q is replaced by ^ /T^- Substituting Eqs. Q and 
© into Eq. ©, we obtain 

= -z/:/")-i{[QQp(")+H.c.] 
-[Q(-V^""'^Q + h.c.] 

-[Q(+V^"+'^Q + H.c.]}. (6) 

Here Q(±) = &^\C)Q, G(±)(£) = /^^ diG(±)(Oe-*^*, 

and Q — + Q'^'. Under the wide-band approxi- 

mation for the detector reservoirs, the spectral function 
&^\C) can be expHcitly carried out as 0: G^^^(/C) = 
1] [x/(l - e"^/'^)]^^_^^^^, where ry = 2TTgLgR, and T is 
the temperature. In this work we will use the unit sys- 
tem oi h = e ~ kB = I- In Eq. ijHJl the terms in {• • • } 
describe the fluctuation effect of the forward and back- 
ward electron tunneling through the detector on the mea- 
sured system. In particular, the Liouvillian operator "£" 
in G^^)(£) characterizes the energy transfer between the 
detector and the measured system, which correlates the 
energy (spontaneous) relaxation of the measured system 
with the inelastic electron tunneling in the detector. At 
high- voltage limit, formally V C, the spectral function 
Ci±)(C) ~ G(±)(0), and Eq. ® reduces to the previous 
result derived by Gurvitz et J, 0,13,01 • 

In the following, we specify the measured system as a 
pair of coupled quantum dots (a solid-state charge qubit) , 
described by the Hamiltonian _ffqu = €a\a){a\ +eb\b){b\ + 
il{\b){a\ + \a){b\). Introduce e (e^ — eb)/2, and set 
(ca + e&)/2 as the reference energy. The qubit eigen- 
energies are obtained as Ei = ve^i-fF = A/2, and 
Eq — — Ve^ -|- fl'^ = —A/2. Correspondingly, the eigen- 
states are |1) = cos ||a) + sin||6) for the excited state, 
and |0) = sin ||a) — cos for the ground state, where 
6 is introduced by cos 6* = 2e/A, and sin 61 = 20/ A. 
The coupling between the qubit and detector is char- 
acterized by H' = QF, where Q ~ T + x\o.){o-\, and 

^ = Em(4A + H-c.)- 

With the knowledge of p''"'\t), one is able to carry out 
the various readout characteristics of the detector. In the 
strong projective measurement regime (e.g. = 0), the 
measurement-induced wavefunction collapse of the qubit 
can be perfectly manifested by the probability distribu- 
tion function P(n,t) = Tr[p'^"^(t)]. Switching on Vl such 
that l/i7 is comparable to or smaller than the measure- 
ment time 0, the qubit state oscillation cannot be read 
out by a series of single shot measurement. In this regime, 
the continuous weak measurement is an alternative ap- 
proach to register the qubit oscillations. In the remained 
part of the paper, we calculate the output current and 
noise spectrum based on Eq. © . 

Straightforwardly, the average current flowing through 
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the detector can be generally expressed as 

lit) = ^nTr[p(")(i)] = iTr[QpQ + H.c.], (7) 



where Q = ^ — Q^^-*. For symmetric qubit (i.e., e = 
or 6 = 7r/2), the stationary current reads 



Is = 9oV + giV 



1 - 



(8) 



Here go = vCT + 5i = ^/(x/2)^ and G(±) = 

i [F(+)(A,y)±F(-)(A,y)], with F(±)(A,y) = (A± 
coth(-^^). We notice that Eq. coincides with the 
result derived in Ref. ^ but differs from that in Ref. 0. 
The former was obtained on the basis of real-time Green's 
function diagram technique, while the latter was resulted 
from the quantum trajectory technique under rotation- 
wave approximation. In addition to the measurement 
current, in the following we detail the studies of output 
noise spectrum in the regime of continuous weak mea- 
surement, where a number of remarkable new features 
will be revealed. 
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FIG. 1: (a) Noise spectrum in the presence (solid curves) and 
absence (dashed curves) of the qubit relaxation, (b) Volt- 
age effect on the noise spectrum, particularly on the peak-to- 
pedestal ratio (inset, where the solid and dashed curves cor- 
respond to the presence and absence of the qubit relaxation). 
The results in (a) and (b) are obtained, respectively, by al- 
tering X (for a fixed voltage V — 2A ) and the voltage V (for 
a fixed X ~ O.IA). Other parameters are Ql = qr = 2.5/A, 
and r = A. 



The noise spectrum can be calculated using the Mac- 
Donald's formula 



S{uj) ^2uj dtsmuut— [{n^{t)) - [Itf] 
Jo 



(9) 



where / is the average current over time, and (n^it)) 
n'^Pin, t). It can be shown that 



dt 



QN{t)Q + ^Qp{t)Q + B..c. 



(10) 



where N{t) = J2n''^P v-)' which can be calculated via 
its equation of motion 



dN , 1 

— = -iCN - - 
dt 2 



Q,QN-NQ^ +i(QpQ + H.c.)(.ll) 



For symmetric qubit, it would be desirable to carry out 
the exphcit result. Denoting S{uj) = So + Si{uj) + S2{i^), 
the result reads 



So = 2/oCoth 



2T 2 



G(+) - - V^coth — 
G(+) 2T 



Si{uj) 



A G(- 



2V G(+) 



(tj2_ A2)2+r2cj: 
G(-) 



^2 + r2- 



(12a) 
r,(12b) 
(12c) 



Here three currents are defined as Iq = {la + Ib)/2, Id = 
h-Ib, and J- /o- ir;x'AG(-)/G(+), with /„ = ry(r + 
x)^V and lb — rfT^V being the detector currents corre- 
sponding to qubit states \a) and |6), respectively. Other 

quantities in Eq. p2|l are introduced as: Yd — ^^G^+\ 

7 = 2|iA and = -A^/Jj;/G^+^ - rix^G^-^A. The 
three noise spectrum components are, respectively, (i) the 
zero-frequency noise So, (ii) the Lorentzian spectral func- 
tion Si{uj) with a peak around the qubit Rabi frequency 
= A, and (iii) 5*2 (w) completely originating from the 
qubit relaxation induced inelastic tunnelling effect in the 
detector. In addition to 5*2 (w), the qubit relaxation also 
manifests its effect in So and 5*1(0;), i.e., giving rise to the 
second term of So and reducing the pre- factor in 6*1(0;) 
from unity. If the qubit relaxation induced inelastic effect 
is neglected, or at the limit of high bias voltage y ^ A, 
Eq. I12II returns to the known result of previous work 
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The measurement-induced relaxation effects of the 
qubit are further shown in Fig. 1. The major effect of 
the qubit relaxation shown in Fig. 1(a) is lowering the 
entire noise spectrum, in qualitative consistence with the 
finding by Gurvitz aZ where an external thermal 
hath is introduced to cause qubit relaxation. However, 
the spontaneous relaxation discussed here does not di- 
minish the telegraph noise peak near zero frequency in 
the incoherent case, which implies the surviving of the 
Zeno effect, in contrast to the major conclusion of Ref. 
0. Also, the transition behavior from the coherent to the 
incoherent regime is different. Figure 1(b) shows the volt- 
age effect that the coherent peak around w = A reduces 
as decreases the measurement voltage. Interestingly, this 
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effect alters the fundamental upper bound limit of 4 for 
the value of the peak-to-pedestal ("signal-to-noise") ra- 
tio, [S{A) — So]/ So, which was found by Korotkov et al. 
at the high voltage limit (see the inset) 0. 



(a) e = 71 12 (b) e = 71 16 




FIG. 2: 3D-plot of the noise spectra for (a) the symmetric 
qubit, and (b) the asymmetric qubit. The adopted parameters 
are c/l = pj? = 2.5/A, x = O.IA, and T = A. 

The voltage effect is further shown in Fig. 2 by the 3D- 
plot of the scaled spectra for different qubit symmetries. 
In contrast to the present result, we notice that in Ref. 10 
no spectral structure was found, i.e., S{uj) — S{oo) — 0, 
in the wide range oi V < lOA for the symmetric qubit 
{6 = tt/2). However, Shnirman et a? showed the existence 
of the coherent peaks at uj = ±A for voltage higher than 
A 0. For asymmetric qubit as shown in Fig. 2(b), the 
coherent peaks at cj = ±A are destroyed and a peak 
around w = is formed. This transition originates from 
the breakdown of the resonant condition, which replaces 
the Rabi oscillation of the qubit by incoherent jumping. 

Finally, simple analysis in limiting cases can provide 
additional insight into the correlation between the detec- 
tor and the qubit. At zero temperature, we obtain a pre- 
factor 1/2 in Si{oj) in the low- voltage regime {V < A). 
This result is in sharp difference from previous conclu- 
sions: In Ref. i such kind of contribution vanishes, while 
in Ref. ^3 it does not exist at all in a much wider range 



of voltage. At the same limit, Eq. H12() also predicts non- 
vanishing 5*2(0;) and non-zero correction to the Schottky 
shot noise 2/o in So- Remarkably, all these contribu- 
tions were in absence from the previous work '9', Td] , and 
the reason was attributed to the complete relaxation of 
the measured qubit to its ground state. Here we under- 
stand our distinct result as follows. Consider the key 
quantity {n'^{t)) — J^n^'^^i''^'^)- Despite the fact that 
the qubit would relax to its ground state under the con- 
cerned limit, the fluctuation of (n^(i)) remains according 
to Eq. lO, since the conditional qubit state p'^^\t) is not 
at all stationary. In particular, the non-zero off-diagonal 
elements of /o*-"-* (i) contain the information of qubit co- 
herence, which gives rise to the peak structure of the 
noise spectrum. Therefore, differing from the previous 
work ^ JLO I and even going beyond the very recent de- 
bate |ll| , we conclude here that at zero temperature and 
even in low bias voltage, the detector output noise spec- 
trum still contains excess components in addition to the 
Schottky shot noise, due to the quantum fluctuations in- 
duced by the coupling of the detector and the qubit. This 
novel feature may deserve further confirmation in future 
work. 

In summery, we have generalized the validity range of 
the quantum measurement theory developed by Gurvitz 
et al. to arbitrary voltage and temperatures. The gen- 
eralized theory properly accounts for the energy transfer 
between the detector and the measured system. Its ap- 
plication to charge qubit measurement reveals a number 
of interesting new readout characteristics associated with 
the new treatment of the correlation between the detec- 
tor and the qubit. 
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